Charged domain walls as quantum strings living on a lattice. 



o 

Q 



CO 

c3 



C 
O 

> 

m 
o 

&\ 

-(— » 
c3 



c 
o 
o 



X 



Henk Eskes*, Osman Yousif Osman, Rob Grimberg, Wim van Saarloos and Jan Zaanen 
Instituut-Lorentz, Leiden University, P.O. Box 9506, NL-2300 RA Leiden, The Netherlands 

(February 1, 2008) 

Recently experimental evidence is accumulating that the the doped holes in the high-T c cuprate 
superconductors form domain walls separating antiferromagnetic domains. These so-called stripes 
are line-like objects and if these persist in the superconducting state, high-T c superconductivity is 
related to a quantum string liquid. In this paper the problem of a single quantum meandering string 
living on a lattice is considered. A kink model is introduced for the string dynamics which allows 
us to analyze lattice commensuration aspects. Building on earlier work by den Nijs and Rommelse, 
this lattice string model can be related both to restricted Solid-on-Solid models, describing the 
worldsheet of the string in Euclidean space time, and to one-dimensional quantum spin chains. 
At zero-temperature a strong tendency towards orientational order is found, and the remaining 
directed string problem can be treated in detail. Quantum delocalized strings are found whose 
long wavelength wandering fluctuation are described by free field theory, and it is argued that the 
fact that the critical phase of delocalized lattice strings corresponds to a free Gaussian theory is 
a very general consequence of the presence of a lattice. In addition, the mapping on the surface 
problem is exploited to show the existence of new types of localized string phases; some of these 
are characterized by a proliferation of kinks, but the kink flavors are condensed so that the long 
wavelength fluctuations of these strings are suppressed. The simplest phase of this kind is equivalent 
to the incompressible (Haldane) phase of the S = 1 spin chain, and corresponds with a bond centered 
string: the average string position is centered on bonds. We also find localized phases of this type 
which take arbitrary orientations relative to the underlying lattice. The possible relevance of these 
lattice strings for the stripes in cuprates is discussed. 



I. INTRODUCTION 



A series of experimental developments has changed the 
perspective on the problem of high T c superconductivity 
drastically. As long as the doping level is not too high, 
electrons bind at temperatures well above T c |lj, and 
the superconducting state appears to be in tight com- 
petition with some collective insulating state j|. There 
exists compelling evidence that this insulating state cor- 
responds with a novel type of electron crystal, character- 
ized by both spin- and charge condensation: the stripe 
phase [§-§[. This phase consists of a regular array of 
charged magnetic domain walls: the holes introduced by 
doping form line- like textures which are at the same time 
anti-phase boundaries, separating antiferromagnetic spin 
domains — see Fig. la. This stripe phase is observed 
in systems where the insulating state is stabilized by Zn 
doping pUq] or by the so-called LTT collective pinning 
potential P,0- 

Inelastic neutron scattering data reveal that strong dy- 
namical stripe correlations persist in the metallic- and su- 
perconducting regimes HJq~H] ■ Although no static stripe 
order is present, the magnetic fluctuations as measured 
by inelastic neutron scattering should reflect stripe corre- 
lations. As was shown very recently, the magnetic mod- 
ulation wavevector of the static stripe phase seems iden- 
tical with that of the dynamical spin fluctuations in the 
metal- and superconductor for various doping levels M. 



In addition, it was argued that the anomalous normal 
state magnetic dynamics can be explained in terms of 
domain wall meandering dynamics M. 

The exciting possibility arises that the zero temper- 
ature superconducting state is at the same time a rel- 
atively mildly fluctuating quantum stripe fluid. Unlike 
the rather featureless diagonal sector of, e.g., A He [ |10[ , 
it can be imagined that the charge- and spin sectors of 
this quantum stripe problem have an interesting internal 
structure. Because charged domain walls are line-like ob- 
jects, the charge sector might be looked at as a quantum 
string liquid |S| Jll| , |l2| . Little is known in general about 
such problems, and theoretical analysis is needed. In or- 
der to address the problem of many interacting strings, 
it is first necessary to find out the physics of a single 
string or charged domain wall in isolation. A string is 
an extended object, carrying a non-trivial collective dy- 
namics — in contrast to particle-like problems, the el- 
ementary constituent of the string liquid poses already 
a serious problem. The physics of quantum strings is 
a rich subject. This is most easily discussed in terms 
of path-integrals. In D+l dimensional Euclidean space 
time, a particle corresponds with a world line, and so the 
quantum string corresponds with a "worldsheet" . The 
statistical physics of membranes is a rich subject, which 
is still under active investigation [Of . 

The debate on the microscopic origin of the stripe in- 
stability is far from closed Q-ffl,p|,p|. Nevertheless, in 




this paper we will attempt to isolate some characteristics 
which might be common to all present proposals for the 
microscopy, to arrive at some general considerations re- 
garding the quantum meandering dynamics. From those 
we will abstract a minimal model for the string dynam- 
ics. The phase diagram of this model can be mapped out 
completely, and turns out to be remarkably rich. 
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FIG. 1. (a) A charged domain wall separating spin do- 
mains of opposite AFM order parameter, (b) Breaking up 
domain walls causes spin frustration, while (c) "kinks" do 
not. (d) Kinks can gain kinetic energy by moving along the 
domain wall, (e) A typical rough wall, (f) An example of a 
directed string. 

These characteristic features are: (i) It is assumed that 
the charge carriers are confined to domain walls. This is 
the major limitation of the present work and it is hoped 
that at least some general characteristics of this strong 
coupling regime survive in the likely less strongly coupled 
regime where nature appears to be. (ii) In addition, we 
assume that domain walls are not broken up, as sketched 
in Fig. 16, as this would lead to strong spin frustration. 
(in) Most importantly, we assume a dominant role of lat- 
tice commensuration on the scale of the lattice constant. 
Configuration space is built from strings which consist of 
"holes" which live on the sites of an underlying lattice. 
An example of such a string configuration is sketched in 
Fig. lc. This automatically implies that the microscopic 
dynamics is that of kinks along the string (Figs. 1 c,d), 
and this leads to major simplifications with regard to the 
long wavelength behavior of the string as a whole. Note 
that there is ample evidence for the importance of lattice 
commensuration: the scaling of the incommensurability 
with hole density x for x < 1/8 0, the special stability 



at x = 1/8 j|, the LTT pinning mechanism ||. (iv) It 
is assumed that the strings do not carry other low lying 
internal degrees of freedom, apart from the shape fluctu- 
ations. Physically this means that localized strings would 
be electronic insulators. The data of Yamada et al. H in- 
dicate that this might well be the case at dopings x < 1/8 
(the linear dependence of the incommensurability on x in- 
dicates an on-domain wall charge commensuration), but 
it is definitely violated at larger dopings where the strings 
should be metallic fll9|-f22||. Work is in progress on fluc- 
tuating metallic strings, where we find indications that 
the collective string dynamics is quite similar to what is 
presented here for insulating strings [|23|. 

Given these requirements, one would like to consider a 
quantum sine Gordon model J24| for the string dynamics, 
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Here z(l) is the transversal displacement at point I on the 
string, and II(Z) its conjugate momentum defined through 
the commution relation [II(Z), z(l')] = i8(l — I'), and c is 
the transversal sound velocity. The first two terms in 
Eq. (1.1) describe a free string, while the last term is 
reponsible for the lattice commensuration effects: every 
time the string is displaced by a lattice constant, the 
potential energy is at a minimum. This model is well 
understood J24J . When the strength of the nonlinear in- 
teraction exceeds a critical value (g > g c ), the interaction 
term is relevant and the string localizes. The excitation 
spectrum develops a gap and it is characterized by well- 
defined kink and anti-kink excitations. When g < g c the 
sine term is irrelevant, and although the dynamics is at 
least initially kink-like on microscopic scales, the string 
behaves as a free string at long wavelength. The latter 
is the most elementary of all quantum strings. It follows 
immediately that the relative transversal displacement of 
two points separated by an arclength I along the string 
diverges as ((z(l) — z(0)) 2 ) ~ InZ j|). The string as a 
whole is therefore delocalized, and this is the simplest 
example of a "critical" string. 

A central result of this paper is that Eq. flLl| ) is, at 
least in principle, not fully representative for the present 
lattice problem. More precisely: starting from a more 
complete microscopic kink dynamics model (section II) 
we find a richer infrared fixed point structure. The phase 
diagram incorporates phases associated with the quan- 
tum sine Gordon model fixed point, but also includes 
additional phases which are intimately connected with 
the effects of the lattice and of the nearest neighbor in- 
teractions between the holes. In section III, we derive 
the path integral representation of our model. It turns 
out that the worldsheet of this string in Euclidean space 
time corresponds with two coupled restricted solid-on- 
solid (RSOS) surfaces |25y, each of which describes the 



motion of the string in either the x or y direction on the 
two dimensional lattice. 

The bare model is invariant under rotations of the 
string in space. As discussed in section IV, we find 
indications for a generic zero-temperature spontaneous 
symmetry breaking: for physical choices of parameters, 
the invariance under symmetry operations of the lattice 
is broken. Even when the string is critical (delocalized 
in space), it acquires a sense of direction. On average, 
the trajectories corresponding with the string configu- 
rations move always forward in one direction while the 
string might delocalize in the other direction, see Fig. 
1/. This involves an order-out-of-disorder phenomenon 
which relatively easy to understand intuitively. Quan- 
tum mechanics effectively enhances the fluctuation di- 
mension by stretching out the the string into a world 
sheet in the time-wise direction, and the enhancement of 
the effective dimension reduces the effect of fluctuations. 
Thermal fluctuations destroythis directedness, but they 
do so more effectively when the string isless quantum 
mechanical. 

This directedness simplifies the remaining problem 
considerably. We will show that the directed string prob- 
lem is equivalent to a well known problem in surface sta- 
tistical physics: its worldsheet is equivalent to a single 
RSOS surface. At the same time, this model is easily 
shown to be equivalent to a generalized XXZ quantum 
spin chain problem. The particular model we study is 
actually equivalent to the 5 = 1 spin chain, which has 
been studied in great detail. The RSOS surface problem 
and the quantum spin chain problem are therefore also 
related to each other. This equivalency was actually at 
the center of the seminal work of den Nijs and Rommelse 
on the hidden order in Haldane spin chains pq ]. From 
our perspective, the introduction of the physically ap- 
pealing quantum string model as an intermediate model 
which connects both with the spin chain and the RSOS 
surfaces, also helps to appreciate the depth of the work 
of den Nijs and Rommelse P6[ . 

The bulk of this paper (sections V-VIII) is devoted to 
an exhaustive treatment of this directed string model. 
Some powcrfull statistical physics notions apply directly 
to the present model, and these allow us to arrive at a 
complete description of the phase diagram of the quan- 
tum string. As was announced already in a short com- 
munication |DJ, this phase diagram is surprisingly rich: 
there are in total ten distinct phases. In the context of 
the quantum spin chain/RSOS surfaces, already six of 
those were previously identified. However, viewing this 
problem from the perspective of the quantum string, it 
becomes natural to consider a larger number of poten- 
tially relevant operators and the other four phases be- 
come obvious. 

Compared to strings described by Eq. (1.1), we find 



energy so that the string as a whole is localized — we use 
"localized" here in the sense that the transversal string 
fluctuations of two widely separated points remain finite, 
((z(£) — z(0)) 2 ) — ► const, as I — > oo. Besides the differ- 
ent directions the purely classical strings can take in the 
lattice, we also find a number of localized strings which 
have a highly non-trivial internal structures: the "dis- 
ordered flat" strings, characterized by a proliferation of 
kinks, but where the kink flavors condense so that the 
string as a whole remains localized. On the other hand, 
the quantum-delocalized (critical) strings are all of the 
free field variety and as we will argue in the final section, 
this might be a very general consequence of the presence 
of a lattice cut-off. 



II. MODEL: THE MEANDERING LATTICE 
STRING 



Whatever one thinks about the microscopy of the 
stripes, in the end any theory will end up considering the 
charged domain walls as a collection of particles bound 
to form a connected trajectory, or such a model will be 
an important ingredient of it. Moreover, these trajecto- 
ries will communicate with the crystal lattice, because 
the electrons from which the strings are built do so as 
well. This fact alone puts some strong constraints on the 
collective dynamics of the charged domain walls. 

Let us consider the string configuration space. On the 
lattice this will appear as a collection of particles living 
on lattice sites, while every particle is connected to two 
other particles via links connecting pairs of sites. The 
precise microscopic identity of these particles is unimpor- 
tant: they might be single holes (filled charged domain 
walls |l4[]l5f l as in the nickelates |27[]), an electron-hole 
pai r (the charge density waves of Nayak and Wilczek 



or a piece of metallic- 
domain wall. All what 



a much richer behavior but this is limited to the regime 
where lattice commensuration dominates over the kinetic 



[12 1, or Zaanen and Olcs 
[20 1 or even superconducting [2 
matters is that these entities have a preferred position 
with regard to the underlying lattice (site ordered |L4j, 
or bond ordered |19| ). Quite generally, curvature will 
cost potential energy and a classical string will therefore 
be straight, oriented along one of the high symmetry di- 
rections of the lattice. Without loss of generality, it can 
be assumed that the lattice is a square lattice while the 
string lies along the (1,0) ('#') direction. Denoting as N y 
the number of lattice sites in the y direction and assum- 
ing periodic boundary conditions, this straight string can 
be positioned in N y ways on the lattice. Obviously, such 
a string will delocalize by local quantum moves: the par- 
ticles tunnel from site to site 0,E3|. Moving the whole 
string one position in the y direction involves an infinity 
of local moves in the thermodynamic limit, and the dif- 
ferent classical strings occupy dynamically disconnected 
regions of Hilbert space. 



This is analogous to what is found in one dimen- 
sional systems with a discrete order parameter |3Cj. In 
the case of, e.g., polyacetylene the order parameter is 
of the Zi kind: the bond order wave can be either 

A-B-A-B or B-A-B-A 

(A single bond, B double bond), while a single transla- 
tion over the lattice constant transforms the first state 
of the staggered order parameter into the second kind 
of state. This is a discrete operation, because the lat- 
tice forces the bond-order to localize on the center of the 
bonds. Such an order parameter structure implies the 
existence of topological defects, which are Ising domain 

walls: A-B-A-B-B-A-B-A ("kink") and 

B-A-B-A- A-B-A-B ( "antikink" ) . When 

they occur in isolated form, these are also genuine build- 
ing blocks for the quantum dynamics, because although 
their energy is finite, it involves an infinity of local moves 
to get rid of them (topological stability) . In the particu- 
lar problem of polyacytelene, these kinks only proliferate 
under doping (charged solitons). Although topological 
quantum numbers are no longer strictly obeyed when the 
density of topological defects is finite, it has been shown 
in a number of cases that they nevertheless remain gen- 
uine ultraviolet quantities as long as they do not overlap 
too strongly [fl|,|||. 

If we consider a (locally) directed piece of string, the 
string is analogous, except that the symmetry is now 
Zn : on the torus, a half-infinity of the string is local- 
ized at y position n y , and the other half can be displaced 
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1. Hence, in total there are 



N y — 1 distinct kink excitations with the topological in- 
variants corresponding with the net displacement of the 
half-string in the y direction. Because the kink operators 
can occur in many flavors, this problem is therefore in 
principle richer than that of one dimensional solids. 

Clearly, kinks with different flavors have to be dynam- 
ically inequivalent. Since there is apparently a reason 
for the particles to form connected trajectories, it should 
be more favorable to create a kink corresponding with a 
small displacement than one corresponding with a large 
jump. Here we will focus on the simplest possibility: 
only kinks occur, corresponding with a displacement of 
one lattice constant in the y-direction. This restriction is 
physically motivated by the fact that the string is thought 
to separate two antiferromagnetically ordered states; so, 
if the displacement of successive holes would be larger 
than one lattice constant, the antiferromagnetic ordering 
would be strongly suppressed — after all, this is the very 
reason that holes tend to line up in stripes. In addition, 
we will specialize on the "neutral" string. It will be as- 
sumed that the string is characterized by a gap in its 
charge and spin excitation spectrum, so that the strings 
with kinks contain the same number of particles as the 
classical reference configurations. The model we will con- 
sider might apply literally to the charge commensurate 
stripes of the nickelate (27fl. In the cuprates, it might be 



better to consider the stripes as one dimensional metals 
or superconductors, characterized by massless internal 
excitations. In these cases, it remains to be demonstrated 
that eventually the transversal string fluctuations decou- 
ple from the internal excitations for the present model to 
be of relevance. 

Given these considerations, we propose the following 
model for quantum lattice strings. The string configu- 
rations are completely specified by the positions of the 
particles ("holes") r; = (xi,yi) on the 2D square lattice. 
Two successive particles I and I + 1 can only be nearest or 
next-nearest neighbors, or |r;+i — r;| = 1 or y/2. We will 
call these connections between successive particles links. 
Two classes of links, those of length 1 and those of length 
y/2, exist. Taking the order of the particles into account 
there are 8 distinct links. The string Hilbert space is 
spanned by all real space configurations satisfying the 
above string constraint. 

We consider local discretized string-tension interac- 
tions between nearest and next-nearest holes in the 
chain,(H = -0H) 

Hex = J2 [*2(l*J+i " x i\ ~ l)*(lM+i - Vi\ - 1) 



Y^ Cij5(\xi+i -Xi-i\ -i)S(\yi +1 -yi-i\ - j)] 



i,3=0 



M Y S ( T l ~ T m) 



(2.1) 



The various local configurations and interaction energies 
are shown in Fig. 0. The last term is an excluded vol- 
ume type interaction — the physically relevant limit is 
M — ► oo, so that holes cannot occupy the same site. 
The interaction K, distinguishes horizontal from diago- 
nal links, and Cij = Cjt is a set of two- link interactions, 
which one can think of as microscopic curvature terms. 
Furthermore, we exclude strings with a physically un- 
realistic extreme curvature by taking Cio — ► oo. Note 
also that configurations which would give an contribu- 
tion Coo to the energy are automatically excluded in the 
limit M — > oo, which we will take throughout this paper. 
There are five local configurations, distinguished by four 
parameters. Therefore we can choose £20 = and the 
string is determined by the parameters /C,£n,£i2, and 
£22, see Fig. §. 
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FIG. 2. (a) The set of local configurations and their clas- 
sical energies, (b) The two allowed hopping processes. We 
take t — t'. 

The second term in 7i is a quantum term which allows 
the particles to hop to nearest neighbor lattice positions. 
However such hopping processes should not violate the 
string constraint. This constraint can be enforced by 
means of a projection operator Pstr(^i+i — rj), that re- 
stricts the motion of hole I to the space of string config- 
urations, 



PstAr) = 6(\r\ - I) + 6(\v\ - V2) 



(2.2) 



The string is quantized by introducing conjugate mo- 
menta n a , [r™,7T£,] = i8i_ m 8 a ^ 1 where a — x or y. A 
term e m7r < acts like a ladder operator, and causes parti- 
cle / to hop a distance n in the x-direction, 



\xi) = \xi +n). 



(2.3) 



Therefore the kinetic energy term becomes, 
Wo = 2r53i^ P (r I+1 -r I )Pf tP (r,-r,_ 1 )coB(7rf). 
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(2.4) 



Note that the particles, even under Hq, keep their order 
and therefore can be labeled by I. Thus the fermion na- 
ture of holes in realistic domain walls plays no role at our 
level of approximation, and quantum statistics becomes 
irrelevant. 

The above model is minimal since it contains only 
nearest neighbor hopping and the simplest string ten- 
sion terms. One natural extension would be to take the 
two hopping amplitudes T in Fig. g& to be different, 
since there is no microscopic reason why they should be 
identical. In the next sections we will discuss the zero 
temperature properties of the above string model. The 
self- avoidance term is a complicated non-local operator. 
However, we will find that, surprisingly, the kinetic en- 
ergy favors oriented walls without loops. Therefore this 
term turns out to be unimportant for the present zero- 
temperature discussion. 



III. RELATION TO RSOS-LIKE SURFACE 
MODELS 

The problem introduced in the previous section can 
be reformulated as the classical problem of a two dimen- 
sional surface (worldsheet) embedded in 2+1 dimensional 
space, using the Suzuki- Trotter mapping. The model can 
be seen as two coupled RSOS (restricted solid-on-solid) 
surfaces. The solid-on-solid models are classical mod- 
els for surface roughening p5[. They describe stacks of 
atoms of integer height in two dimensions, with an inter- 
action between adjacent stacks depending on the height 



differences. With this construction overhangs are ex- 
cluded. In the restricted SOS models these height differ- 
ences are limited to be smaller or equal to some integer n. 
In the present case, the two RSOS models parametrize 
the motions of the world sheet in the spatial x- and y 
directions, respectively, while the (strong) couplings be- 
tween the two takes care of the integrity of the world 
sheet as a whole. 

In the Suzuki- Trotter [|33| or Feynman path integral 
picture one writes the finite temperature partition func- 
tion as an infinite product over infinitesimal imaginary 
time slices. In this limit the commutators between the 
various terms in the Hamiltonian vanish like l/n 2 , where 
n is the number of Trotter slices, and the partition func- 
tion can be written as, 



Z = lim Tr 



(3.1) 



To show the relation with RSOS models, we will cast 
the transfer matrices T in the form of a two-dimensional 
classical effective Hamiltonian. This implies writing the 
matrix elements of the T-matrix between configurations 
{]";} in terms of an effective classical energy depending on 
the worldsheet positions {r;^.}, where k is the imaginary 
time index running from 1 ton with periodic boundary 
conditions. Schematically, 
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Since Tici is diagonal in the real-space string basis, it is 
already in the required form, 
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For Tiq a few more steps are needed, 
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The expression in the last line is of course only valid for 
states in which the aj's in successive time slices differ by 
at most one unit. Combining these two energy contribu- 
tions we arrive at the following classical problem, 
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This classical world sheet is constrained to \xi t k+i — 
Xl,k\ < 1 an d \yi,k+i — yi,k\ < 1, and the interactions 
are anisotropic. The above classical model can be viewed 
as two coupled two-dimensional RSOS surfaces, xi t k and 
yi t k- The x coordinate of hole I at the time slice k is 
now identified as the height of an RSOS column posi- 
tioned at (I, k) in the square lattice. In a similar way 
the y coordinates define a second RSOS surface, coupled 
strongly to the first by the above classical interactions. 
Since the steps Ax can at most be equal to 1, the RSOS 
sheets are restricted to height differences 0, ±1 between 
neighboring columns. 

The classical model as defined above is not unique. 
While the above mapping allows us to exploit the con- 
nection to other models most efficiently, for the numeri- 
cal Monte-Carlo calculations a different decomposition is 
used, which allows for a more efficient approach to the 
time continuum limit. We write, 



Z = T re Hc '+' H Q 
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(3.6) 



In the above formula I is the identity operator, in our 
case a complete set of string configurations. We have 
chosen to split the T-matrix into a contribution from 
even and odd sites, or A and B sublattice (checkerboard 
decomposition) , 
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and a similar expression for the odd sites. 7i(2l) is the 
Hamiltonian of the even string element 21, equal to Eq. 
(2T) or (|]J) without the sum over string links. L is the 
number of links in the chain. Because of the sublattice 
decomposition, Ta is a simple product of local T-matrices 
and Z becomes, 

n 1/2 

z=lim E [Hl« r 'WA l H r '}i) 

{r,, k y ik }k=ll = l 

({r,}' fe |t 2 J +1 < fe |{r,} fe+1 ). (3.8) 



Each timeslice is split in two subslices, r and r'. The no- 
tation {ri} k denotes the set of positions r; at the given 
time slice with index k. Note that the t matrices are in- 
dependent of / and k , and these indices only label the 
position of the t matrix in the 2D world sheet. The local 
t matrices, t A and ts, depend only on three positions. 
For instance, 



({r,} fe |^' fc |{r,}' fe ) 
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(3.5) witn ^e restriction r 2 i-i,k = *2i-i,k anci r 2i+i,k 



r 2(+i fe- Since each link has 8 different orientations, the 
local t matrix connects in general 8 x 8 = 64 possibilities. 
However most of the t matrix elements are zero, and it 
decomposes into subblocks, of which the biggest one is 
3x3. The states which are connected via the local t 
matrix, or the Hamiltonian, are listed in Fig. 0, and the 
matrix elements are given in the appendix. The above 
decomposition will be used in the Monte-Carlo simula- 
tions described in Section IV. 
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FIG. 3. The four subblocks of the local t-matrix. The 
other equivalent, symmetry related blocks are obtained by 
7r/2 rotations and reflections in the x or y axis. 



IV. DIRECTEDNESS AS SPONTANEOUS 
SYMMETRY BREAKING 



We are not aware of any similarity of the statistical 
physics problem of the previous section with any exist- 
ing model. RSOS problems are well understood, but it 
should be realized that in the present model the two 
RSOS problems are strongly coupled, defining a novel 
dynamical problem. When we studied this problem with 
quantum Monte-Carlo, we found a generic zero temper- 
ature symmetry breaking: although the string can be 
quantum delocalized, it picks spontaneously a direction 
in space. This symmetry breaking happens always in the 
part of parameter space which is of physical relevance. 



Our understanding of this phenomenon is a nice illustra- 
tion of an order-out-of-disorder mechanism which at first 
sight appears paradoxical, but which is easy to under- 
stand qualitatively. At first sight, one might expect that 
the quantum fluctuations (kinetic energy) would tend to 
disorder a string, i.e., to decrease the tendency for the 
string to be directed. That the opposite effect happens, 
is due to the fact that in the presence of quantum flucua- 
tions, the string traces out a worldsheet in the imaginary 
time direction. The larger the kinetic term, or the larger 
the temperature, the further the worldsheet stretches 
out in the timewise direction. At zero temperature, the 
worldsheet becomes infinite in this direction as well. The 
statistical physics of a string is then equivalent to that of 
a fluctuating sheet in three dimensions. Now, it is well 
known from studies of classical interfaces [Q that while 
a one-dimensional classical interface in two dimensions 
does not stay directed due to the strong fluctuations, 
for a two-dimensional sheet the entropic fluctuations are 
so small that interfaces can stay macroscopically flat in 
the presence of a lattice |3^,|36| - - for this reason, the 
roughening transition in a three-dimensional Ising model 
is properly described by (i.e., is in the same universal- 
ity class as) a Solid-on-Solid model in which overhangs 
are neglected pqpfl|. In other words, even if microscopic 
configurations with overhangs are allowed, a classical in- 
terface on a lattice in three dimensions can stay macro- 
scopically flat or "directed" . In agreement with this pic- 
ture, we will find that the directedness order becomes 
more robust when the string quantum-fluctuations be- 
come more severe, i.e., when the wordsheet stretches out 
further into the timewise direction. 

It is convenient to use open spatial boundary condi- 
tions. In this case, the string configuration space of the 
model defined in section II consists of all surfaces, con- 
necting to the boundaries in the space-directions while 
periodic along the imaginary time direction. On a time 
slice, the trajectories can be categorized according to the 
global property of directedness: when the string connects 
two opposite sides of the lattice in the spatial directions, 
this trajectory can be called "directed" and otherwise it 
is undirected. An alternative, more restrictive definition 
of directedness is: every continuous string configuration 
s can be written as a parametrized curve in two dimen- 
sions, [x(t),y(t)], where t could for instance be the dis- 
crete label of the successive particles along the string. 
When the string configuration can be parametrized by 
a single- valued function x(y) or y(x), we call the string 
configuration directed (see Fig. p]/) . The quantum string 
vacuum is a linear superposition of many string config- 
urations. When all configurations in the vacuum cor- 
respond with single valued functions x{y) or y(x), the 
string vacuum has a directedness order parameter. We 
will give an explicit measure of the directedness in our 
simulations below. 

It should be immediately clear that directedness or- 



der is rather fragile. It cannot exist at any finite tem- 
perature. When temperature is finite, the width of the 
worldsheet in the imaginary time direction becomes finite 
as well, and the long wavelength fluctuations of the string 
becomes a ID statistical problem, subjected to string 
boundary conditions which are inconsequential for the 
argument which follows. Consider first the classical limit: 
at zero temperature the string would be straight, running 
along (say) a (1,0) direction. A local 'corner' configura- 
tion of the type shown in Fig. Ha would be an excitation 
with energy L\\ (alternatively, one could consider two 
kinks). Clearly, a single corner suffices to destroy the 
directedness of the classical ground state. At any finite 
temperature, the probability of the occurrence of at least 
one corner is finite: P = A r cxp(— f3Lu). Hence, direct- 
edness order cannot exist at finite temperatures, for the 
same reasons that long range order is destroyed at any 
non-zero temperature in one dimension. 

At zero temperature, the internal dynamics of the 
worldsheet is that of a two-dimensional sheet which is 
infinite in both directions and since the lattice commen- 
suration has rendered the problem to be of a discrete 
symmetry, long range order can exist for finite kinetic 
energy. We will discuss the symmetry breakings of a sin- 
gle RSOS surface in great detail later. For the present 
discussion it suffices to know that such a single RSOS sur- 
face can be fully ordered, as well as (partly) disordered. 
Because of the strong coupling, it would a priori appear 
questionable to discuss the dynamics of the full model of 
Section II in terms of the dynamics of the two separate 
RSOS subproblems. However, in the context of direct- 
edness it is quite convenient to do so. When both the 
x and y RSOS problems would be fully disordered, it is 
easy to see that the string vacuum would be undirected. 
This is illustrated in Fig. ^a: two kinks moving the string 
from a (1,0) to a (0,1) direction in the lattice correspond 
with one kink which can move freely in the horizontal 
part of the string, and one kink which can move freely 
in the vertical part of the string. On the other hand, 
when both RSOS problems are ordered, the string is also 
ordered. For instance, the (1,0) string can be thought of 
as a combination of an RSOS surface which always steps 
upwards in the x direction, and one which is horizontal 
in the y direction (Fig. |4p). 

A third possibility is that one of the RSOS subprob- 
lems is ordered, while the other is disordered. Dismissing 
crumpled phases (like condensates of the L\\ type cor- 
ners), the only possibility remaining is that one of the 
RSOS problems steps up always, while the other is dis- 
ordered, as illustrated in Fig. He. This results in a dis- 
ordered directed string vacuum: the string steps always 
forward in, say, the x direction while it freely fluctuates 
in the y direction. We find that this symmetry breaking 
happens always as long as the string does not crumple. 
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FIG. 4. (a) An undirected string with two kinks propagat- 
ing along different directions. Note that the bend blocks the 
propagation of kinks, (b) a (1,0) string and the corresponding 
two (coupled) RSOS surfaces along the x and the y directions 
respectively. The numbers correspond to x (y) position of 
hole I at imaginary time t. (c) a disordered directed string 
and the corresponding ordered and disordered RSOS surfaces 



In the latter case half of the string is rotated around 
any of the sites I G 2..L — 1. This means that for 
instance the position of all holes m > I are replaced 
by (x m ,y m ) ->• {xi,yi) + (y m - yi,-{x m - x t )). Such 
operations turn out to be very efficient — completely 
wrapped high temperature strings unwrap in just a cou- 
ple of Monte Carlo steps at low temperature. 

We discovered this directedness symmetry breaking 
during the simulations. To quantify the directedness we 
measured the fraction of string length in the string vac- 
uum corresponding with single valued configurations. At 
zero temperature, the ground state wave function of the 
string is, 



l*o) = Yl a ({xi,yi})\{xi,yi}) 



(4.1) 



{xi,yi} 



where every state in string configuration space (\{xi, yi})) 
corresponds with a trajectory [x(t),y(t)). Consider first 
the case of a continuous string. For every configuration, 
the total string arclength is given by 



L{{x u yi})tot 



ds 



\/dx 2 + dy 2 



(4.2) 



To perform Monte-Carlo calculations on the string 
model, we used the checkerboard decomposition as de- 
scribed in Sec. III. Note that the string has no minus 
sign problem — the sign of T is easily transformed away 
by an appropriate choice of the operators. For £ 10 — ► oo, 
there are 8x5=40 different local configurations, given by 
the positions rj_i, r; and rj + i. However most of the 
40 x 40 local T matrix elements are zero, and the largest 
subblock is 3 x 3. These blocks were listed in Fig. ||. 
The corresponding matrix elements are listed in the ap- 
pendix. For the Monte-Carlo program to produce sensi- 
ble results it is crucial to have operations that add and 
remove bends easily. We added global mirror and 7r/2 
rotation operations illustrated in Fig. a. 
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FIG. 5. Two additional Monte Carlo operations used for 
the simulations of the general string. The 90 degree rotation 
(a) around position 1 in this example turns a non-directed 
string into a directed one. The mirror operation (b) is im- 
portant to quench defects in diagonal strings (mirror plane 
indicated by the dashed line). 
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FIG. 6. Illustration of the way we measure the direct- 
edness of a string in the continuum case (a) and on the lat- 
tice (b). The heavy solid parts of the string indicate the 
parts where the projection of the string onto the x axis is sin- 
gle-valued, and for which the indicator function g y (x) equals 
1. 

To quantify the directedness, we introduce, as illus- 
trated in Fig. o, an indicator function g y (x) which equals 
1 when the string is single-valued when projected onto 
the a>axis, and zero otherwise, and analogously a func- 
tion g x (y) which equals 1 when the curve is single- valued 
when projected onto the y-axis, and zero otherwise. The 
total directed length in the x and y directions are then 
defined as 



L ({xi,yi})dir,x = dx g y {x) 




L({x h yi})dir, y = / dy g x (y)\jl + ( -j- 



(4.3) 



The measure of directedness in the string vacuum is then 
defined as the larger of N% ir (Q) and N% ir (Q), where 



^(0) 



{xi,yi} 



\a\{x hyi })\ MxuVl})d ^ 



L{xi,yi)tot 



(4.4) 



where n = x,y. On the lattice, our measure of direct- 
edness is the immediate analogue of this definition, ex- 
cept that we just count the number of directed bonds, 
irrespective of whether they are oriented diagonally or 
horizontally. 

By thermal averaging, the above definition of direct- 
edness density is immediately extended to finite temper- 
ature, 



NL(T) 



£< 



-/3(E„-E ) 



NIM) 



(4.5) 



where N dir (n) is the directedness density of a string ex- 
citation with energy E n . Eq. (4.5) can be straight- 
forwardly calculated using quantum Monte-Carlo. A 
Monte-Carlo snapshot defines a stack of coupled string 
configurations along the imaginary time direction (the 
Trotter direction). We calculate Nf ir for every Trotter 
slice by calculating the fraction of the string length in this 
configuration which is single valued in the ^-direction. 
This is given by the number of bonds which steps for- 
ward in the x-direction devided by the total number of 
bonds in the string. We then average this quantity over 
the string wold sheet (Trotter direction) and then over 
the Monte-Carlo measurements. The same is done for 
N v dlr {n). The larger of N^ r (n) and N v dir {n) is then the 
density of directedness at the given temperature. 

That the zero-temperature properties of the directed 
and non-directed string are equivalent for Cn > can be 
understood by the following argument. The 7r/2 bends 
in strings block the propagation of links along the chain. 
Close to the bend itself the particles in the chain cannot 
move as freely as in the rest of the chain. This effect is 
shown in Fig. 0. 



(K) 





<KKK) CK>00 

FIG. 7. Illustration of the fact that a bend blocks the 

propagation of links along the string. Note holes 1 and 2 

cannot move. 



In space-time the 7r/2 bend is like a staight rod in 
time. Therefore the presence of such kinks increase the 
kinetic energy. For the argument it makes little difference 
whether the bend consists of a single 7r/2 corner or two 
7r/4 corners. This confirms that it is the kinetic energy 
which keeps the strings oriented along one particular di- 
rection. In terms of a directedness order parameter this 
result implies that such a quantity is always finite, except 
when C\\ < or when the hopping term vanishes (it is 
easy to see that in the classical case, T = 0, in many 
regions of parameter space the problem becomes that of 
a self-avoiding walk on a lattice in the limit T — ► 0) . For 
the two equivalent RSOS surfaces this means that one 
of the two spontaneously orders and becomes diagonal, 
while the other RSOS sheet shows any of the ten quan- 
tum rough or classically ordered phases discussed later. 

In Fig. g we show results of typical Monte-Carlo cal- 
culations for the density of directedness as a function of 
temperature Ndi r (T). We have considered four points 
in the parameters space; as will be discussed later, 
these points are representative for phases with interesting 
quantum fluctuations, and serve to clarify our conclusion. 
In Fig. ||(a) the triangles (dashed line) is the result for 
the density of directedness at the point where all the clas- 
sical curvature energies are zero, i.e., corresponding to 
the pure quantum string. The crosses (dotted line) and 
the filled squares (dashed-dotted line) are the results for 
points where /C = 1.8 and K. = 4.0 respectively, and the 
rest of the classical curvature energies are zero. In term 
of the phase diagram for the directed string problem of 
Fig. pThn section VI and Table III in section VII, the first 
point corresponds to a Gaussian string (pure quantum) 
and the other two correspond to flat strings. The point 
K, — 1.8 lies just inside the flat string phase II where sig- 
nificant quantum fluctuations are still present, while the 
point K, — 4.0 lies deep inside the flat phase. The fourth 
curve in Fig. |^(a), given by the full line, is the result of 
a Monte-Carlo calculation for a classical string (T = 0) 
where only flat segments and ir/2 corners are allowed (no 
diagonal segments). This same classical result is shown 
again in Fig. ||(b) together with the result of the direct- 
edness density for a point in the midlle of the Gaussian 
(XY) phase [K = 0.5, £ 2 i = -0.25, C 22 = -1-0, C u = 
corresponding to D — and J = —0.5]. 

We will first discuss the result for the classical case, 
a classical string (T = 0) with Cn, the energy of the 
7r/2 corner, equal to one. Because the string is of finite 
length, the infinite temperature limit of Ndir (T) is not 
zero but rather a small but nonzero value [B7| (~ 0.03 
for a domain wall of length 50). Ndi r (T) is already close 
to this value for all temperatures of order Cn and larger. 
For an infinitely long domain wall Ndi r {T) drops very 
fast to zero with increasing temperature. At the other 
limit, for low T where T <C Cn, Ndi r (T) grows very fast 
to 1. 
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FIG. 8. A Monte-Carlo result for the directedness density 
Ndir{T) at 4 points, (a) The XY point (Triangles) where all 
curvature energie are zero. Two points in the flat phase, with 
K, — 1.8 (crosses) and K, — 4.0 (filled squares) the rest of the 
curvature energies are zero, (b) Inset: A point in the middle 
of the gaussian phase with parameters tC = 0.5, £21 = —0.25, 
£22 = —1.0, £11 = (open circles). The full line in both 
figures is the result for a classical string where only flat bonds 
and -| corners are present with £n = 1. 

Again because the string is of finite length, it becomes 
directed already at a finite temperature: For all temper- 
atures such that Lexp(— (3C\i) < 1 the string configura- 
tions in our simulations are typically completely directed. 
An infinitely long classical string becomes directed only 
at T = 0, of course, since at any nonzero temperature al- 
ways some corners will occur in a sufficiently long string. 

For the quantum string, all the curves look strikingly 
similar to the classical one. When the temperature is 
very much higher than the kinetic term, T 3> T, all 
curves merge together and the classical limit is reached. 
At low T, where r«T, Ndi r (T) again grows very fast 
to 1, as in the classical case; it reaches this value at a 
finite temperature for the finite length string. This is 
even true for the purely quantum string at the XY point, 
where all classical microscopic curvature energies are zero 
(see the dashed line in Fig. H(a)). We can understand this 
in terms of an effective corner or bend energy C that is 



produced by the quantum fluctuations. As in the clas- 
sical case the probability for the occurance of a bend is 
proportional to ~ exp(— /?£). At zero temperature no 
bend is present and the string becomes directed. A fi- 
nite length string effectively becomes directed already at 
a temperature such that Lexp(— (3C) < 1. At intermedi- 
ate temperatures, where the temperature is of the order 
of the kinetic term, things are more difficult and it is far 
from obvious what is going on. Especially in this region, 
all the various classical curvature energies may play a 
role, and the interplay of these on the directedness is un- 
clear. Nevertheless, as is clear from the data of Fig. 0(a), 
this region connects the high and low temperature lim- 
its smoothly. Moreover, by comparing the results for the 
three quantum strings in this figure it is also clear that 
when the string is more quantum mechanical Ndi r (T) is 
higher. 

Our general conclusion, based also on Monte-Carlo 
studies of the behavior in other phases which are briefly 
summarized below, is that apart from some extreme clas- 
sical limits, the general lattice string model at zero tem- 
peratures is a directed string. The qualitative picture of 
7r/2 turns blocking the propagation of kinks appears to 
be a natural explanation for these Monte-Carlo findings. 
The phase diagram of the general string model intro- 
duced in section II will essentially be the same (apart 
from special limits) as the corresponding phase diagram 
of the simplified directed string model. In the remain- 
ing sections of this paper we will therefore focus on the 
phases and phase transitions of the directed string. 

For completeness, we end this section with a brief qual- 
itative description of our observations concerning spon- 
taneous directedness at low but finite temperatures in 
regions of the phase diagram where the directed string 
has other type of ordering than that already discussed. 
All the results apply to L\\ — 0, and we refer to Table 
III in section VII for a quick introduction to the various 



phases of the directed string problem and for the num- 
bering (I-X) of the varous phases. 

- The entire zero temperature phase diagram of the di- 
rected string is reproduced. 

- Phase I is very stable with respect to bends. With 
"stable" we mean that finite strings do not change their 
appearance when increasing the temperature from zero 
to a moderately small temperature, of the order of 0.1 T . 

- Deep in the horizontal phase II (large positive /C) 
quantum fluctuations are strongly suppressed, and at the 
same time the string becomes susceptible to tt/2 corners. 
On the other hand, when we approach from phase II the 
boundaries with phases IV and V, the fluctuations in- 
crease and the string stiffens, Fig. H This is in agreement 
with the picture sketched before that quantum fluctua- 
tions orient the string. 

- Deep inside phase III the string changes constantly 
between horizontal zigzags and vertical zigzags. A tt/2 
turn costs no extra energy. Again close to phase V quan- 
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turn fluctuations have the effect of removing bends. 

— The Haldane phase V and the rough phase IV are very 
robust, and a considerable fraction of 7r/2 bends occurs 
only at relatively high temperatures of the order of 0.2 
T. 

— In the slanted phase VII high temperatures are needed 
before down diagonal links come in. On the other hand 
horizontal links are easily replaced by vertical ones. This 
only increases the energy very slightly, but the entropy 
gain is considerable. A typical low temperature string is 
shown in Fig. O in section VII. To zeroth order the hor- 
izontal and vertical links can be thought of as spinless 
fermions moving coherently along the string. In the di- 
lute limit these links have only a weak interaction. The 
order of the links is conserved, and at zero temperature 
the ground state has only horizontal links. However our 
simulations indicate that for a small range of negative 
Ln values a diagonal string with alternating horizontal 
and vertical links is favoured. It is again the kinetic en- 
ergy of the horizontal and vertical links that keeps the 
string oriented in the (1,1) direction. 




FIG. 9. A typical low temperature string in the slanted 
parameter region VII. 



position of a single "guider point" . As a guider point we 
can take the position r of any one of the particles, which, 
together with the relative coordinates given by the links, 
fixes the position of the entire string. Since the guider 
represents just a single degree of freedom, and since the 
thermodynamic behavior of a chain is determined by the 
link interactions, the guider coordinates will be irrelevant 
for the behavior of the chain. Apart from this guider de- 
gree of freedom the directed string problem reduces to a 
one-dimensional quantum problem with three flavors. 



From Eq. (3.5) one directly deduces the Hamiltonian 



of the string directed along x, 



n eff = z2 



/./,■ 



K 



S(\yi+i,k -yi,k\ - 1) 



c 



12 



— S(\yi+i,k -yi-i,k\ - 1) 
n 

h s (\yi+i,k - yi-i,k\ - 2) 

n 

T 
+ ln(- )S( \y, tk - y itk+1 - 1 

n 



(5.1) 



It is clear that the directedness simplifies the model con- 
siderable. The directed version can not self-intersect, and 
the excluded volume constraint is satisfied automatically. 
Furthermore the £n-type of configurations are not al- 
lowed, thus the directed model is specified by three pa- 
rameters and the temperature (T = 1). Because of the 
preceding considerations, Eq. (5.1) corresponds with a 
1+1 D problem, which is actually equivalent to a general 
quantum spin-1 chain. 

We identify the spin with the string height difference 
Vi+i — Vii which can be either 0, 1 or —1, see Fig. nol 



DIRECTED STRINGS AND THE SPIN-1 
CHAIN 



Quite generally, the string problem does not simply re- 
duce to that of the internal dynamics of the worldsheet, 
because of the requirement that the worldsheet has to 
be embedded in D+l dimensional space. However, in 
the presence of directedness order and in the absence 
of particle number fluctuations |2l| the string boundary 
conditions are trivially fulfilled and the string problem 
is equivalent to that of a single "world sheet " in 1+1 
dimensions. Assume the string to be directed along the 
x direction. Since the string steps always forward in this 
direction, the number of particles in the string has to be 
equal to the number of lattice sites in the ^-direction, 
and every directed string configuration will connect the 
boundaries in this direction. The string is still free to 
move along the y direction. Instead of labeling the po- 
sitions in the 2D plane the string is completely specified 
by the list of links, for which there are only 3 possibil- 
ities (in the (1,1), (1,0), or (1,-1) direction), and the 
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FIG. 10. The relation between spin-1 and directed strings, 
Sf = yi+i -yi- 

These link dynamical variables specifying the string 
can be directly identified with the m s = 0,±1 variables 
of the spins living on the sites of the spin chain. Defining 
the latter using hard core bosons b\ riB , the spin opera- 
tors for the S — 1 case become S z = b\b\ — &Li&-ij S + = 
y/2(b[bo+b b-i) and by comparing the action of the spin- 
and string operators on their respective Hilbert spaces 
one arrives at operator identities |26|. A quantum hop 
from y to y + 1 increases the height difference on the left 
of I by one, and decreases it by one on the right, as is 
easily seen by inspecting the two hopping terms in Fig. 
0. Therefore fl2l, 
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S" = yi+\- vi 
St_ x Sf = 2P str (yi - m-i)Pstr(yi+i - yi) 

The following identities, for S = 1, 



,±»*i 



(5.2) 



<5(|y i+1 - yi |-i) = (Sf) 2 , 
S(\vi+i w-i| - 1) = (5f ) 2 + (5f_ x ) 2 - 2(Sf 5f_ x ) 2 , 

*(|l/i+i ~ W-il - 2) = -Sf SU[1 + S?Sf +1 ] , (5.3) 

are easily checked. The directed string problem can now 
be reformulated in spin language as, 



H 



Spill 



Ei 



c 



22 



(/C + 2£ 12 )(Sfr + -^5fSf_ 1 



Co 



l (=f-2C ia )(S?S?_ 1 ) 2 



7" 



(5.4) 



Following the spin-1 literature we define the following 
parameters, 



D = 1C + 2£ 12 , 

J = £12/2, 

-E = £22/2 — 2£i2 



(5.5) 



The -E term is new. It is a quartic Ising term, leading 
to extra phases and phase transitions. For the special 
choice E = (T = 1), the above Hamiltonian reduces to 
the familiar XXZ model with on-site anisotropy, 



Hxxz = ]T [D(Sf ) 2 + JSf ^ 



"2 ^ ^-1 + S i ^i-i) 



(5.6) 



The zero temperature phase diagram of the above spin- 
1 model has been discussed in detail in the literature. 
p8|- [4l| , p6[ In the next section we will briefly review the 
six phases found for this model, from a string perspective. 
Then we will show that a nonzero E parameter leads to 
the appearance of four extra phases in section VII. 

Den Nijs and Rommelse |26| discuss a direct mapping 
between the spin chain and the RSOS surface. We stress 
that this mapping in fact involves two steps: First the 
RSOS model is mapped on a string problem, using the T 
matrix. Then the spins are identified as shown above. 
Thus the quantum string is a natural intermediate of 
the two other models. Den Nijs and Rommelse make 
use of the freedom in the choice of the T matrix to de- 
fine a mapping which is slightly different from ours, since 
they introduce a transfer matrix along a diagonal, while 
we introduce one along the rr-direction. As a result, in 
their case there are only interactions between next near- 
est neighbors along the (1,1) direction, while our choice 



allows for interactions between next nearest neighbors 
along the x-direction. Therefore, our RSOS model dif- 
fers slightly from theirs. 

The RSOS representation is more transparent than the 
quantum model. The spin-1 phases and the nature of 
the phase transitions all have a natural interpretation 
in space-time. For instance the Haldane phase, or AKLT 
wavefunction, with its mysterious hidden string order pa- 
rameter is identified as a "disordered flat" RSOS surface 
pa with a simple local order parameter. The height rep- 
resentation, dual to the spins, gives a similar local order 
parameter for the quantum string. 



VI. THE PHASES (E=0) 

In this section the general string Hamiltonian will be 
simplified by leaving out the quartic Ising term [E = 
in Eqs. ( p 7 ^ ) and (|5.5D). Our string problem is now equal 
to the spin-1 XXZ model. The zero-temperature phase 
diagram of the string problem is surprisingly rich, and 
even for the case E = there are 6 phases and a large 
variety of phase transitions. These phases can be classi- 
fied in three groups: classical strings localized in space, 
quantum rough strings of the free variety, and partly de- 
localized phases of which the disordered flat phase is a 
remarkable example. In this section we will briefly review 
the six phases as discussed in the literature on the spin-1 
XXZ problem (5.6). The problem will be addressed from 
the quantum string perspective. For more details we re- 
fer to Ref. [£6|. In the next section we will show that 
with a finite E > four additional phases are stabilized. 

The phase diagram of the quantum string is shown in 
Fig. [LI 1 as a function of D and J. We ha ve used the 
XXZ model parameters, defined in Eq. (5.5), such that 



the phase diagram can be compared directly with the 
spin-1 literature |3q-pd[| and in particular with Fig. 13 
of Ref. |2(| . We will introduce below the various order 
parameters that have been introduced in this reference 
to distinguish the six phases in this phase diagram. The 
relation between the more general [E 7^ 0) string and 
spin phases will be clarified in the next section. 

There is first of all a horizontal and a diagonal string 
phase. In the diagonal phase I no quantum fluctuations 
are allowed, since, as shown in Fig. || a diagonal string 
does not couple to other states by Hq . This phase is sta- 
bilized by a large and negative £22, so that since E = 
also J = £12/2 = £22/8 is large and negative. A suitable 
variable introduced to define order parameters, following 
Ref. (p6|, is the Ising spin variable 07 = (— l) Vl , which 
identifies whether a given height is in an even or odd 
layer. This underlying spin model can have "ferromag- 
netic" or antiferromagnetic" order, and so we introduce 
the corresponding order parameters [|42J 
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Pstr = (<n(yi+i -vi)) 



(6.1) 



order parameter p s t a g, which is zero in phase IV. In this 
phase horizontal links are virtual and occur in pairs. As 




FIG. 11. The phases and phase transitions of the directed 
quantum string as a function of the on-site anisotropy D 
and the Ising interaction J of the corresponding spin-1 XXZ 
model. The parameter E is set to zero. 



Here the brackets denote the ground state expectation 
value as well as an average over string members I. In 
(3.1), we have also included the order parameter p str dis- 



cussed below. In the horizontal phase II one particular 
height is favoured, thus the order parameter p is nonzero 
here. This phase is stabilized by a large positive /C, which 
suppresses diagonal links. However TLq causes virtual 
transitions from two horizontal links into two diagonal 
ones, see Fig. 0. On the 2D worldsheet these fluctuations 
show up as local terraces that do not overlap and thus 
do not destroy the long-range order. In both phases the 
elementary excitations are gapped. 

Upon lowering K, the terraces grow and at some point 
they will form a percolated network — the string has 
become disordered in both space and imaginary time. 
Via the well-known Kosterlitz-Thouless roughening tran- 
sition |35|] phase IV is entered for J < 0. This phase 
belongs to the well known XY universality class, char- 
acterized by algebraic correlation functions and gapless 
meandering excitations — capillary waves in fluid inter- 
face language. The roughness, however, is extremely 
"soft" and the height difference diverges only logarith- 
micly, ((yi — J/m) 2 ) ~ In |/ — m|. The transition from the 
Gaussian phase which is rough and on average oriented 
horizontally to the "frozen" diagonal phase is a "quasi 
first-order" KDP transition |28]. 

For large negative K, diagonal links are favored over 
horizontal ones. There is a transition to a second rough 
phase (phase VI) . It is distinguished from the first by the 



we will discuss later in section VII , for large negative K 
the model can therefore be reduced to an effective spin 
1/2 problem. 

For negative K, and positive J {—L\ij2 — £22/8) the 
string becomes a (physically unlikely) zigzag with alter- 
nating up and down diagonal pieces. Excitations to pairs 
of horizontal links are gapped. Again p s tag = (( — 1) Zcr ?) 
serves as an order parameter. Upon increasing K. the 
islands formed by pairs of horizontal links start to over- 
lap and there is an Ising transition into the Haldane or 
disordered flat (DOF) phase. 

The point J = 1,D = belongs to the gapped 
DOF phase, in agreement with Haldane's educated guess 
[]43| , |44 1 that integer spin chains are gapped at the Heisen- 
berg antiferromagnetic point. In this "disordered hori- 
zontal" string phase the prototype wavefunction, equal 
to the AKLT valence bond state p5[ , has every up diag- 
onal link followed by a down link, with a random number 
of horizontal links in between. The height yi takes just 
two values, say and 1. The local order parameter is, 
p s tr, defined in (pT|). This order parameter measures the 
correlation between the next step direction and whether 
one is in a layer of even or odd height. When p s t r = 1, 
the string just steps up and down between two layers, but 
the steps can occur at arbitrary positions. Note that the 
height is a global quantity in spin language, i.e., it is the 
accumulated sum over spins, yi = J2 m =o ^m- Because of 
this the above order parameter becomes non-local when 
rewritten in terms of the "string" of spins. Therefore, it 
is often called the string order parameter. We will also 
use this name, but stress that the "string of spins" to 
which this name refers should not be confused with the 
general strings which are the basis of our model, and that 
the other order parameters are nonlocal as well in terms 
of the original spins S. 

This phase diagram can be rationalized by writing 
the RSOS problem as the product of a 6-vertex model 
and the 2D Ising model of s spins on the 6-vertex 
lattice, as discussed in detail by Den Nijs and Rom- 
melse. (26|. The horizontal, diagonal, zigzag but also 
the second rough phase VI all correspond to Ising or- 
der: p — (ai) is nonzero in the horizontal phase II, while 
Pstag = ((— 1) 07) is nonzero in the diagonal phase I, the 
zigzag phase III and the rough phase VI. The six-vertex 
part is defined on the crossing points of steps on the sur- 
face — see Fig. O. This is a (sometimes highly) diluted 
set of points. The Ising degree of freedom disorderes on 
the transition between the phases III and V, and between 
IV and VI, while the six-vertex part remains unchanged. 
Therefore these transitions are Ising like. Transitions I 
-»• IV, I -> VI , IV -> V and III -> VI are related to the 
six-vertex part becoming critical, and these KDP and 
KT transitions are known from the quantum spin-1/2 
chain. The transition II to IV is related to the famous 
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surface-roughening transition, of the Kosterlitz-Thouless 
type. |35|,[36| The subtle transition between phase II and 
V, is coined "preroughening transition" by den Nijs. It 
separates two gapped phases. At the transition the gap 
closes and the system is Gaussian, with varying expo- 
nents along the transition line, pq-ffll 
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FIG. 12. Vertices (thick dots) on the space, imaginary 
time string worldsheet. The numbers correspond to the 
heights yi t h. Arrows are drawn when the heights of neigh- 
bors differ. When four arrows occur at a crossing point this 
is called a vertex. 

Almost all the phases can be distinguished by the 
above order parameters p, p s tag and p s t r , except that 
these do not discriminate between the diagonal phase I 
and the rough phase VI. These two phases can be identi- 
fied by also introducing an order parameter which detects 
the presence of an average slope, p s i op e = (yi+i — Hi)- In 
Table Q we list the various phases for E = and the order 
parameters. 

As we shall see in the next section, in the general case 
E 7^ it is more convenient to introduce slightly differ- 
ent spin variables to identify all the ten different phases 
that occur then. The choice of Ref. [gfj] discussed here 
is somewhat more convenient for understanding the uni- 
versality classes of the various phase transitions. 



VII. THE FULL PHASE DIAGRAM - 
PHASE-BOUNDARY ESTIMATES 

As mentioned above the quartic Ising term with pref- 
actor E generalizes the XXZ Hamiltonian and leads to 
extra phases. We will motivate that four extra phases are 
to be expected, and we will show that they are stabilized 
by a positive E parameter. The most disordered phase 
is still the Gaussian phase. 



Using a similar decomposition as above, we can deter- 
mine how many different phases to expect for a general 
spin-1 chain with z-axis anisotropy and nearest-neighbor 
interactions |4q] . Think of the spin 1 as consisting of two 
spins 5, see Table O. The first is a z —\ when the spin 1 
has S z = and a z =| when S z = ±1, similar to the Ising 
degree of freedom defined above. This spin thus indicates 
the presence or absence of a step. The second spin i s, is 
defined as s z = S z /2 when S z = ±1 and is absent when 
S z = 0. This is related to the diluted vertex network 
discussed by den Nijs and Rommelse, in that if there is a 
step, the z-component of s indicates wether this step is 
up or down. The spins s can have ferromagnetic (F) or 
antifcrromagnetic (AF) order, or they can be disordered 
(D). For a the two ferromagnetic cases correspond to 
different physical situations, and we have to distinguish 
ferromagnetic J, (Fl), a horizontal string, from ferromag- 
netic | (F2). When a has Fl order, s becomes irrelevant 
(or better — there are disconnected finite terraces of s 
spins with short-range correlations). Therefore one ex- 
pects 10 phases, depending on the order of the two spin 
species: 1 Fl phase, 3 F2 phases, 3 c-disordered phases, 
and 3 cr-antiferromagnetically ordered phases. These are 
listed in Table III . An example of a phase diagram in a 



case in which all ten phases are present is show in Fig. 
|l3| , which corresponds to the case E — 5. The detailed of 
how this phase diagram was obtained will be discussed 
below. 




FIG. 13. The phases and phase transitions of the quan- 
tum string for E — 5 as a function of t. On the axis are the 
on-site anisotropy D and the Ising interaction J. 

There are four new phases, VII to X, compared to the 
phase diagram discussed in the previous section. All four 
are stabilized by a positive E parameter in Eq. (p]j). 
Three phases , VIII to X, result from an antifcrromag- 
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netic order of the a spin. This corresponds to alternating 
horizontal and diagonal string links (see Table III). The 
diagonal links can be either all up (FM, phase VIII), al- 
ternatingly up and down (AFM, phase X) or disordered 
(phase IX). In phase VII the a spin is disordered, while 
the s spin is in the FM phase. This is a diagonal wall di- 
luted with horizontal links. These links coherently move 
up and down along the wall, lowering the kinetic energy. 
The wall can take any average angle between — 7r/4 and 
7r/4, and this angle is fixed by the value of the param- 
eters. We will call this the "slanted" phase. In terms 
of the decomposition into an Ising spin model and a six- 
vertex model of den Nijs and Rommelse it is easy to see 
that the horizontal links change the orientation of the 
Ising spin and act like a Bloch wall. The Ising spin is 
therefore disordered. The six- vertex term is irrelevant 
for the existence of the slanted phase — in the case of a 
single horizontal link, i.e., on the boundary between the 
slanted and diagonal string phase, there are no vertices. 
A large part of the phase boundaries can be estimated 
exactly, almost exactly, or to a fair approximation. Let 
us focus first on the classical phases. The diagonal, hor- 
izontal and zigzag phases have the following energies in 
the classical approximation in which there are no fluctu- 
ations, as is easily verified, 

Ei = L{JC + C22) = L(D + J + E), 
E u » , 
E jn ki LK = L(D - J + E) . (7.1) 

where L is the length of the chain. The first-order tran- 
sitions will therefore occur close to the lines K = —£22 
(D = —J — E) between phases I and II, £22 = ( J = 0) 
between phases I and III, and K, = (D = J — E) be- 
tween phases II and III. These transitions become exact 
in the classical or large-spin limit. 

The transition between phase I and VII, the diagonal 
and slanted phase, can be found exactly. The transition 
is of the Pokrovsky-Talapov, or conventional ID metal- 
insulator type (see, for instance, Ref. Q). The hori- 
zontal link can be seen as a hard-core particle or a spin- 
less fermion, with the parameters determining an effec- 
tive chemical potential. For a critical chemical potential 
equal to the bottom of the band of the hard-core parti- 
cle the band will start to fill up. The transition occurs 
when the diagonal string becomes unstable with respect 
to a diagonal string with one horizontal link added. This 
single link delocalizes along the string with a momentum 
k and a kinetic energy 2Tcos(fc). The minimal energy is 
(L - 1)/C + (L - 2)£ 22 + 2£i2 - 2T, and the transition 
occurs when K, = 2(£i2 — £22 — T) or, with T = 1, when 



I to VII transition: D = -2(J + E + 1) 



(7.2) 



The transition between phase III and V will occur 
when horizontal link pairs unbind in the zigzag back- 
ground. A rough estimate, neglecting fluctuations, is 



obtained by comparing the energy of a single horizon- 
tal link with that of a perfect zigzag. In the same way 
as above we estimate the phase boundary to be close to 
D = 2(J — E — 1). In the same way the transition from 
phase II to V or IV is determined by the energy of a 
single diagonal step in a horizontal wall, which becomes 
favourable when D = 2. This last estimate turns out 
to be very crude, in that it largely underestimates the 
stability of the flat phase. 

For large negative /C the horizontal links are strongly 
suppressed, and the string can be mapped perturbatively 
on a spin 1/2 chain. Identify S z = 1 (diagonal upward) 
with s z =1 and S z = — 1 (diagonal downward) with 
s z =|. Via a virtual (0,0) spin pair (two horizontal links) 
the spins can still fluctuate, (1,-1) — > (0,0) — ► (-1,1). One 
finds, using second-order perturbation theory in T/JC, 

H eff (D - -00) = (4J + J± )^ S f S f +1 



-J± 



~5> 



2 

2T 2 



3± 



\2D + ZE\ 



(7.3) 



Here we subtracted an irrelevant constant term. This 
has the form of the well studied spin 1/2 Heisenberg 
chain with Ising anisotropy. Transitions occur when 
(4 J + j±) = ±j±, or when J = (III to VI) and 
J = -1/\2D + 3E\ (I to VI) (setting T = 1). 

The above estimates seem to suggest that the line 
J = is special. Our numerical results show that it 
describes accurately the transition between III and VI, 
but also the transition between IV and V. This agrees 
with the arguments given by Den Nijs and Rommelse 
p6| that the Kosterlitz-Thouless transition between IV 
and V should occur precisely at the J = line. 

The slanted phase consists predominantly of up diago- 
nal and horizontal links. Neglecting down diagonals alto- 
gether, which turns out to be a good approximation, one 
can again map the string or spin-1 chain on an effective 
spin 1/2 system. Now the relevant degree of freedom is 
the a Ising degree of freedom. Because a =\ (a diagonal 
link) is not symmetrically equivalent to a =1 (a hori- 
zontal link) the spins will feel an effective magnetic field, 
which regulates the density of horizontal links. Rewriting 
Eq. Q gives, 

tuff = Dj^tf + \) 

1 

+(J + £)]>>f + i)(af +1 + i) 
1 

+ r E( ,j / +(T r+i + CT r^ + + i0 - (7-4) 
1 

and, after rescaling and putting T = 1, 
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We// = h J2 of + A J2 °1<>1+1 

I I 



(7.5) 



with the field h — (D + J + E)/2 and Ising coupling 
A = ( J + E)/2. On the line h = the number of up di- 
agonal links equals the number of horizontal links. The 
average tilt angle is thus 22.5° in this approximation. The 
phase diagram of the spin 1/2 chain in the h-A plane was 
discussed by Johnson and McCoy pTy , For h — there 
are three phases. The ferromagnet corresponds to phase 
I, the antiferromagnet with phase VIII, and the gapless 
disordered phase translates to the slanted string phase 
VII. Increasing the field h in the AFM phase will cause 
a transition to the gapless phase with a finite magneti- 
zation. In the approximation that down diagonals are 
neglected, it follows from the results of Johnson and Mc- 
Coy @ that the point A = 1, h = or J = 2- E is the 
point with the most negative value of J where phase VIII 
is stable. For E = (as well as for small values of E) 
this occurs in the positive J side of the phase diagram, 
meaning that phase VIII to X will in fact not be stable: 
for positive values of J, down steps in the original model 
proliferate. To have a phase diagram with all 10 phases 
present we choose E = 5. 




FIG. 14. The various phase transitions, obtained from 
semiclassical estimates, exact arguments and perturbative 
mappings to spin 1/2. 

In Fig. Hj the various phase-boundary estimates given 
above are summarized. The topology of the main part of 
the phase-diagram has now become clear. In the centre of 
the figure for E — 5 the Johnson-McCoy phase diagram 
is inserted. The estimates suggest that at least phase VII 



and VIII are stabilized by taking E — 5. The dotted line 
through phase VIII is the line where the effective field h 
is zero, and the number of diagonal links is (nearly) equal 
to the number of horizontal ones. 

We finally argue that the slanted phase exists in some 
region of the phase diagram for any E > 0. To see this, 
first consider the case E = 0. Along the line D = 0, our 
model (with E = 0) corresponds to the Heisenberg model 
with Ising anisotropy, and the point J = — 1 corresponds 
to the isotropic ferromagnetic Heisenberg point. Along 
the line D = 0, the transition from the "ferromagnetic" 
(diagonal in string language) phase I to the XY phase IV 
therefore occurs at J = — 1. Now, for E = 0, the exact 
location of the line along which phase I becomes unsta- 
ble to the sla nted phase VII, given by D — — 2(J + 1) 
according to ( |7.2| ), goes exactly through the ferromag- 
netic Heisenberg point at J = — 1 as well. The results of 
Fig. 13 of den Nijs and Rommelse |2(| indicate that this 
line then touches the phase boundary between phases I 
and IV right at this point in such a way that for E = 
no slanted phase occurs. If we assume that both phase 
boundaries shift linear in E for E nonzero and small, it is 
clear that the slanted phase must stabilize in some region 
near the point D — 0, J = — 1 for one sign of E, while 
for the other sign the phase must be absent. Physically, 
it is clear that the stabilization of the slanted phase will 
occur for positive values of E, E > 0. 



VIII. NUMERICAL 



To fill in the details of the phase diagram we have 
performed exact diagonalization and finite-temperature 
quantum-Monte Carlo calculations. Ground state prop- 
erties of strings up to 15 holes (spin chains of length 14) 
were obtained using the Lanczos diagonalization method. 
For the Monte-Carlo we used the checkerboard decom- 



position, briefly explained in Sec. |I|. The Monte Carlo 
has the disadvantage that an extra limit to zero tem- 
perature has to be taken, a regime where the updating 
slows down considerably and where it is difficult to judge 
the accuracy. For determining the phase boundaries of 
the directed string we mainly used the Lanczos results 
for the equivalent spin model. On the other hand the 
Monte Carlo space-time worldsheets provide a transpar- 
ent physical insight into the phases, phase transitions, 
order parameters, etcetera. Moreover, the Monte Carlo 
method allows, of course, to treat bigger systems. 

We are in the fortunate situation that the order param- 
eters of the various phases and the universality classes of 
the transitions are known. This offers a variety of ap- 
proaches to determine the critical lines — one can mon- 
itor the finite-size behavior of the order parameter, cor- 
relation functions or the energy level spacings. Typi- 
cally we applied two independent methods to the various 
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transitions. Our aim is to map out the entire phase dia- 
gram with an accuracy of roughly the line thickness in the 
phase diagram. For very accurate estimates other meth- 
ods, notably the density-matrix renormalization group 



treatment of White |48|, are more appropriate. 

An elegant and powerful method is the phenomeno- 
logical renormalization group approach pioneered by 
Nightingale. W9 50 In this approach one consideres an 
infinite strip with a width L, as a finite-size approxima- 
tion to the 2D classical system. At the critical temper- 
ature of the infinite system one expects, from finite size 
scaling, that the correlation length along the strip scales 
like the width of the strip {L\ or L2), 



1^2 



(8.1) 



The infinite strip is solved by diagonalizing the (finite) T 
matrix. The correlation length can be calculated from, 



£ = l/ln(Ai/Ao), 



(8.2) 



where Ao and Ai are the largest and second largest eigen- 
values of the T matrix. 

A finite ID quantum chain can be viewed as a strip 
infinitely long in the imaginary time direction (zero tem- 
perature). In the time continuum limit, writing T — 
exp( rg), the equation corresponding to Eqs. (8.1) and 
©is, 

Li[£i(£i) - E {L X )\ = L 2 [E X {L 2 ) - E (L 2 )}, (8.3) 

for two different string lengths L\, L2 — ► 00 and for 
parameter values at criticality. Here Eq and E\ are 
the ground state and first excit ed e nergy of a quantum 
Hamiltonian H. According to (3.3), a phase transition 
line can be located by studying th e en ergy gap as a func- 
tion of L, while moni torin g when 
tical example of Eq 



3j) is obeyed. A prac- 
3) is shown in Fig. [15]. Successive 
curves for lengths L and L + 2 show two crossing points. 
Extrapolating these crossing points to infinite L gives two 
estimates of the preroughening transition from phase II 
to V. 

The above scaling holds when time and space scale in 
the same way, or (E±(q) — Eq) ~ q z , with a dynamical 
exponent z = 1. This can be checked independently, giv- 
ing a self-consistent justification of the use of Eq. (3.3). 
In a similar spirit one can determine the critical scaling 
of correlation functions of an operator O by monitor- 
ing the lowest energy state with a nonzero overlap with 
O|0). This has been used extensively J38|-ffl[ to study 
the phases and exponents of the Gaussian phases of the 
spin-1 chain. We refer to these articles for more details. 

Another method used to determine second-order phase 
transitions is the Binder parameter [ pi] , which we de- 
fine as (3 — (m 4 )/(m 2 ) 2 )/2, where m is the relevant or- 
der parameter. This quantity tends to 1 in the ordered 



phase, where (m 4 ) w (to) 4 , and approaches in the disor- 
dered phase, where m has a Gaussian distribution around 
(to) =0. In a renormalization group sence the shape of 
the order parameter distribution function becomes inde- 
pendent of the size at criticality. The various curves for 
different sizes of the Binder parameter versus model pa- 
rameters or T should therefore cross at the critical point. 
For instance for the Ising transition between phase III 
and V we take to = (— l) Vl ~ L , the Ising degree of freedom 
introduced before. 
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FIG. 15. Estimate of the preroughening transition be- 
tween phase II and V. The plot shows L[E%{L) — Eq(L)] for 
various lengths L, as a function of the parameter D, with 
J — 0.8 and E — 0. The two crossing points between the 
successive curves form an upper and lower estimate of the 
transition. 

The Kosterlitz-Thouless transition between for in- 
stance IV and V is a subtle one due to the infinite order of 
the transition and the exponential vanishing of the gap. 
Previous studies show a large uncertainty in the posi- 
tion of the transition line. In the entire Gaussian rough 
string phase the system is critical. Height correlations di- 
verge very weakly like G(r) = ((y r — yo) 2 ) oc Cln(r). At 
the KT transition point the prefactor takes the universal 
p6[ value C = 2 /n 2 . We found that this relation is very 
useful to determine the KT transition line — see Fig. nn. 
This relation gives surprisingly good results even for very 
small distances, and is consistent with a KT transition 
at J = 0, as discussed by den Nijs and Rommels e. |26[ 

The phase diagram for E — 5 is shown in Fig. |l3j. For 
this value of E the slanted phase VII is very pronounced. 
The phases VIII to X occur in a small region in the mid- 
dle of the diagram around the line of equal probability 
of horizontal and diagonal links, in the approximation 
that the Hamiltonian can be mapped onto the spin 1/2 
problem (7.5). The rough phases occur at small negative 
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values of J. Phases I, II, III, V, VIII and X are gapped. 
1.2 



(b) 




FIG. 16. Estimate of the Kosterlitz-Thouless transition 
between phase IV and V. At the transition the slope between 
C( T ') — ((Vr ~ J/o) 2 ) and ln(r) approaches 2/ty 2 (the dotted 
lines) . 

The character of the various phases becomes clear by 
looking at the corresponding Monte-Carlo snapshots in 

Fig. m 

In the flat phase one particular height yi dominates, 
and the quantum fluctuations do not percolate in space- 
time. Decreasing K. these quantum fluctuation islands 
grow and when they overlap the string enters phase IV. 
The system is Gaussian rough in both space and time. 
Note the very weak logarithmic meandering — e.g. for a 
string of length 10, Fig. |l^ shows that the mean square 
height fluctuations are only of order 1 near the KT tran- 
sition. Despite this, determining the KT transition from 
G(r) works surprisingly well for the small systems calcu- 
lated. 
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FIG. 17. Monte-Carlo snapshots of (a) the Haldane phase 
(Phase V, D = 0, J = 1,E = 0), (b) the flat phase (Phase 
II, D = 2, J = 0,.E = 0), (c) the rough phase (Phase IV, 
D — 0, J — —0.5, E — 0) and (d) the slanted phase (Phase 
VII, D — —0.75, J = —5,E — 5). Black to white means in- 
creasing height, except for the slanted phase in (d), where 
black denotes a horizontal link and white an up diagonal link. 

The fascinating order parameter of the Haldane phase 
or disordered flat phase becomes transparent when look- 
ing at the worldsheet, using the height representation 
instead of the spin-1 language. Globally the surface is 
limited to two heights only, and is therefore macroscopi- 
cally flat. However in both the time and space direction 
there is a disordered array of up and down steps, with the 
restriction that every step up is followed by a step down, 
when the order is perfect. On flat pieces, however, there 
are local fluctuations (with consecutive up steps or down 
steps) decreasing the value of the order parameter. These 
islands will grow when J — > from the positive side, and 
when they overlap the string becomes rough. Note that 
the rigidity (flatness) is clear when viewing the overall 
structure of the worldsheet in 1+1D. On the other hand, 
from a single time slice one might be tempted to conclude 
that the string is rough. 

The transition from phase I to VII is of the Prokovsky- 
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Talapov type. Such transitions are often discussed in the 
context of commensurate-incommensurate transition of a 
monolayer of atoms on a substrate with a different lattice 
parameter. In the "floating solid" phase such a system 
consists of a set of parallel domain walls with entropic me- 
andering. [^H The similarity with phase VII, illustrated 
in Fig. [Tij is clear. The entropic meandering in the 2D 
classical case is now to be interpreted as the quantum 
motion of hard-core particles (horizontal links) along the 
string. 



IX. DISCUSSION AND CONCLUSIONS 

Motivated by stripes, we have introduced a lattice 
string model for quantum domain walls and mapped out 
its full phase diagram. We find a generic zero temper- 
ature symmetry breaking: the string acquires a direc- 
tion in all cases. The main reason is that bends in the 
string prohibit the quantum transport, or, vice versa, 
the quantum motion of kinks straightens out the string 
(the "garden hose" effect of Nayak and Wilczek @). 
We arrive at the counter-intuitive conclusion that for in- 
creasing kink quantum disorder the orientational pref- 
erence of the string grows. The directed string problem 
which remains appears to be related to a well understood 
surface statistical physics (RSOS) model and simultane- 
ously to a S = 1 XXZ quantum spin chain with single 
site anisotropy. Motivated by the string interpretation, 
we found a number of phases described by this class of 
models which were previously not identified. 

Physically, the phases fall in three main categories: 
classical (flat worldsheet), Gaussian (rough worldsheet) 
and 'disordered flat' phases. The phases are further dis- 
tinguished by the direction they take in the embedding 
space. Besides the flat strings in the horizontal and di- 
agonal directions, we find that the disordered flat phases 
show here a rich behavior. Apart from the known phase 
with horizontal direction, which is associated with the 
incompressible phase of the spin model, we identified a 
new category of disordered flat phases which take, de- 
pending on parameters, arbitrary directions in space (the 
"slanted" phases). 

Although this does not apply to the localized strings, 
we suspect that a strong universality principle might ap- 
ply to the delocalized strings: At least away from the 
phase boundaries to the localized phases, the underlying 
lattice renders the delocalized strings to be described by 
free field theory. The reason is simple: regardless of the 
terms that one adds to the lattice scale action, the prob- 
lem remains of the XXZ kind and the massless phases 
fall into the 1+1D 0(2) universality class. For instance, 
one can add other kink-kink interactions, etcetera, and 
these can be all described by products of S z operators. 
Although these operators determine the nature of the lo- 
calized phases, they turn into irrelevant operators in the 



massless phases. The kinetic sector is more subtle. For 
instance, one would like to release the constraint that 
kinks only occur with 'height flavor' ±1. This means in 
surface language that one partially lifts the restricted- 
ness of the RSOS model, or in spin language that one 
increases the total spin: e.g., S = 2 means in string lan- 
guage that kinks occur describing height differences of ±2 
as well. Although increasing the magnitude of spin has 
an influence on the localized phases, it does not change 
the fact that the massless phase away from the phase 
boundaries is still obeying XY universality. A point of 
caution is that the 'holes' in principle could change their 
order when larger excursions are allowed. However, these 
'exchange loops' are strictly local and therefore irrelevant 
for the long wavelength behavior as long as the string is 
internally an insulator. These could represent more of a 
problem for strings which are internally superconductors 
or metals. 

We also stress that it follows from the arguments of den 
Nijs and Rommelse |2(| that the occurrence of a gapped 
Haldane type phase for strings is not a peculiar feature 
of the spin 1 representation, but a general consequence 
of the existence of further neighbor interactions between 
the holes in strings. 

Do our findings bear any relevance to the stripes in 
cuprates? At the very least, they do bring up some in- 
teresting questions: 

(a) Is the stripe solidification in for instance the LTT 
cuprates |p[ in first instance driven by a single string ef- 
fect or by a collective transition of the string liquid? In 
the end it has to be the latter, since a single string cannot 
undergo phase transitions at finite temperatures. How- 
ever, it can be well imagined that the effect of the LTT- 
pinning potential is to stabilize (1,0) directed stripes over 
(1,1) stripes. In the language of this paper, this amounts 
to an increase of the parameter K which could move the 
stripe from the Gaussian phase into the horizontal flat 
phase. At zero temperature, this would turn individual 
stripes in straight rods which are obviously much eas- 
ier to order than meandering strings. At finite temper- 
atures, this could increase the single string persistence 
length substantially, so that stripe-stripe interactions be- 
come more effective in stabilizing a stripe solid at finite 
temperatures M] . Further work is needed to establish if 
these single stripe transitions are of relevance. 

(b) Do the 'disordered flat' string phases exist? The sim- 
plest disordered flat phase is the horizontal one (phase V) 
corresponding with the Haldane phase of the 5=1 spin 
chain. In string language, this is nothing else than a lo- 
calized string along the (1,0) direction in the lattice which 
is however not site-centered (as phase II) but, on average, 
bond centered. Bond centered stripes show up in the nu- 
merical study of the t-J model by White and Scalapino 
]19[ , which shows that the ground state of this model at 
finite dopings is a stripe phase. A main difference with 
the mean- field stripes is that these t-J stripes are bond 
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centered. In first instance, one could be tempted to think 
that this has a truly microscopic reason: charges in t-J 
prefer to live on links. However, it could also be due to 
a collective string effect — it could be "our" phase V. 
This can be easily established by measuring the appro- 
priate (string) correlators. Is it so that on equal times the 
charges live on sites while the kinks take care of delocal- 
izing the stripes over two lattice rows, or is it so that on 
all times the charges are living on the links? This is obvi- 
ously an important question in the light of recent works 
relating the bond centering via Hubbard-ladder physics 
to superconductivity J2^]. We also notice that there are 
experimental indications for bond-centering in the nick- 
elates |53|] where disordered flatness could possibly also 
play a role. 

(c) If well developed stripes exist in the superconductors 
and/or metals, these have to occur in the form of a quan- 
tum disordered stripe phase, or a 'quantum string liquid'. 
What is learnt in this regard from the present study of a 
single string? A prerequisite for the existence of a quan- 
tum string liquid is that a single string is delocalized. If 
our conjecture that a single critical string is described by 
free field theory turns out to be correct, this amounts to a 
considerable simplification. In Euclidean space time, the 
single free string worldsheet is like a Gaussian membrane 
and a system of strings becomes a system of interact- 
ing Gaussian membranes, embedded in 2+1 dimensions. 
This in turn is like a classical incommensurate system in 
3D which, although barely studied, appears as a tractable 
problem. For instance, it is known that the 3D incom- 
mensurate solid melts at a finite temperature in all cases 
[ p4J . For the quantum case this means that the quantum- 
melting transition will occur at some finite value of the 
coupling constant, which in turn depends on the single 
string quantum fluctuation as well as string-string inter- 
action effects. Investigations addressing this many string 
problem are in progress, profiting from the simple fluc- 
tuation behavior of a single string. 
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APPENDIX. 

In this appendix we give the i-matrix elements of the 
blocks defined in Fig. 0. The local i-matrix at position I 
and Trotter slice k is defined as, 



(rj_i, fc , rj,fc, r i+ i, fe | exp -Hi \r'i_ lk , r{ k , r' l+1 k , ) 



(9.1) 



The matrix elements depend on the positions of three 
members of the string, I — 1, I and I + 1. The po- 
sitions of Z — 1 and I + 1 are required to be identical 



Tl-l,k 



l l-l,k l l—l,k 



in the two Trotter subslices in- 




volved, due to the checkerboard decomposition, but the 
position of member I can be different, leading to off- 
diagonal matrix elements. The matrix elements of the 
t-matrix are easily found by first diagonalizing Tit and 
expanding the basis vectors in terms of the eigenvectors. 
Block A: This block contains three configurations, see 
Fig. g. We use the same order for the states as in the fig- 
ure. Note that only half of the energy K, of the diagonal 
link between I — 1 and I and the link between I and I + 1 
should be contributed to I. The Hamiltonian. 



(9.2) 



is easily diagonalized. The eigenvalues are /C, E+ and 
£L. The i-matrix is, 

in t\2 £13 

t = \ tyi t 2 2 tyi 
^13 t\2 Tn 

tn = ~e K/n + N 2 + e E +' n + N 2 _e E ~' n , 
*12 = Nla+e E +/ n + N 2 _a-e E ~ /n , 
ti 3 - --e K ' n + N 2 + e E +' n + N 2 _e E -/ n , 
t 22 = Nla\e E +l n +N 2 _e E -l n , 

± 2 2 

E±-IC 1 

a± = — ;= — , N± = 



T 



2 + a 2 + 



(9.3) 



Here n is the number of Trotter slices and /C, £12, £22, 
£11 and T are the string model parameters. 

Block B contains two configurations, each with one 
horizontal and one diagonal link. Repeating the above 
procedure one finds, 

/ e D cosh(T/n) e D sinh(T/n) \ 
\ e D sinh(T/n) e D cosh(T/n) ) ' 



D 



— 4- ^1H 

2n n 



(9.4) 
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Block C contains a single configuration of two diagonal 
links, and the energy and i-matrix therefore contain £22, 



, ,K £22 x 

t = exp( 1 ), 

n n 



(9.5) 



Block D consists of a square corner between one hori- 
zontal and one vertical link, and C\\ is involved, 



£ = exp( ), 

n 



(9.6) 
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TABLE I. Order parameters that distinguish between the 
six different phases in the phase diagram for E — 0. A + en- 
try in the table indicates that the particular order parameter 



is nonzero. 



Phase p p s tag Pstr Pslope 

I + + 

II + 

III + + 

IV 

V + 

VI + 



TABLE II. Spin 1 S seen as a combination of two spins 
1/2, a and s. 



s 



1 







a 
s 



1 



TABLE III. A schematic representation of the different 
phases. Also shown is the long-range order of the two spins 
1/2, s and a as defined in the text. F = ferromagnetic, Fl = 
up-spin ferromagnetic, F2 = down-spin ferromagnetic, AF = 
antiferromagnetic, D = disordered 



Phase 



I 
II 

III 

IV 

V 

VI 



Fl 

F2 

Fl 

D 
D 

Fl 



String 



AF 

D 
AF 

D 



Spin 1 



++++++++ 
000 00000 



+0-+0+0- 
-+0-0 04 



VII 



D 



0+0++0+0 



VIII 



AF 



0+0+0+0+0 



IX AF D 

X AF AF 



0+0+0-0+0 
0+0-0+0-0 
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